c 
^ 



(N 



Enhanced Saturation Coverages in Adsorption-Desorption 

Processes 

cK 
o\ 

0\ • Paul R. Van Tassel^, Pascal Viot^, Gilles Tarjus^, 

Jeremy J. Ramsden^ and Julian Talbot^ 
^Department of Chemical Engineering and Materials Science, 
Wayne State University, Detroit, MI 48202 
0\ [ "^ Laboratoire de Physique Theorique des Liquides 

Universite Pierre et Marie Curie 
4, place Jussieu 75252 Paris Cedex 05 France 
Q . ^Department of Biophysical Chemistry, Biozentrum, CH-4056 Basle, Switzerland 

^ ' ^Department of Chemistry and Biochemistry, 

Duquesne University, Pittsburgh, PA 15282-1503 

-)— > 

a. 

I . Abstract 

O, 

O , Many experimental studies of protein deposition on solid surfaces involve 

alternating adsorption/desorption steps. In this paper, we investigate the ef- 
feet of a desorption step (separating two adsorption steps) on the kinetics, 
Q>^ . the adsorb ed-layer structure, and the saturation density. Our theoretical ap- 

CN , proach involves a density expansion of the pair distribution function and an 

vQ I application of an interpolation formula to estimate the saturation density as 

^P '_ a function of the density at which the desorption process commences, pi, and 

the density of the depleted configuration, p2. The theory predicts an enhance- 
ment of the saturation density compared with that of a simple, uninterrupted 
^ RSA process and a maximum in the saturation density when p2 = oPi- The 

theoretical results are in qualitative and in semi-quantitative agreement with 
'^ ' the results of numerical simulations. 

o 
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I. INTRODUCTION 

The adsorption of proteins and colloids from solution onto solid surfaces often involves 
strong departure from equilibrium behavior ^. Many experimental studies have attempted 
to quantify the kinetics [||-|^ of the adsorption process, as well as the structure of the 
adsorbed configurations 0,0|. Protein adsorption has also been the subject of a number of 
theoretical studies. In particular, the Random Sequential Adsorption (RSA) model is able 
to account for a number of experimentally observed properties of protein adsorption, such 
as irreversibility and non-linear surface blockage [^?,?]. Quantitative expressions are now 
available for the kinetics of RSA of spherical [^^-|T^ and non-spherical particles [ p!^P 



on 



planar interfaces and these models have been shown to provide an excellent description of the 
adsorption of proteins P,^,P^ , colloids [T^JT^ and other particles. More recently, the RSA 



models have been extended to allow for the possibility of post-adsorption conformational 
change JT^?]. 



It has been observed experimentally that a fraction of the proteins in an adsorbed layer 

may be removed by rinsing with a buffer solution or by eluting with a surfactant solution 

18| , p!9| ,^j3| j20|j2T[1 . The influence of a desorption step on subsequent adsorption is not well 



understood. In particular, one would like to know how the structure of the adsorbed layer 
is affected by the desorption process and how this influences the kinetics when adsorption 
resumes. One anticipates, due to the non-equilibrium nature of the adsorption, that the 
intermediate desorption step may result in behavior that differs non-trivially from an unin- 
terrupted process. It is the goal of this work to investigate this by considering a multistep 
RSA process that begins with the irreversible adsorption of spherical particles onto a planar 
surface up to a density pi. Adsorbed particles are then removed randomly until the density 
falls to p2- We then initiate a second irreversible adsorption process beginning from the 
depleted configuration and investigate how the final saturation density depends on pi and 
P2. The memory effect of the RSA process is expected to result in a final density that is 
different from that of an uninterrupted RSA process. 

In the theoretical part of the study, we use a distribution function approach to determine 
the influence of the multistep adsorption process on the structure of the adsorbed particle 
configuration. A density expansion of the pair distribution function leads to an approximate 
description of the structure up to the maximum density in the initial adsorption step. We 
then show that if particles desorb randomly from the adsorbed configuration, the radial dis- 
tribution function is unchanged. We use this result to estimate the available surface function 
of the configurations generated by the re-adsorption process and develop approximate inter- 
polant formulas that allow us to estimate the final saturation density. The theory is then 
compared with numerical simulations of the process in one and two dimensions. Finally, 
we use these expressions to determine the extent of desorption needed to optimize the time 
required to reach a given adsorbed density. 

II. THEORY 

We consider a three-step adsorption/desorption/readsorption process. The initial phase 
consists of the irreversible addition of spherical particles of diameter a on an initially empty 
planar interface to a density pi < 0.547 (in this work we will use p to denote a dimensionless 



density, p = ja'^N/A where A^ is the number of particles adsorbed in an area A). During 
this phase the particles are placed according to the usual RSA rules i.e., a trial position 
is selected from a uniform random distribution. If the trial particle does not overlap with 
any previously placed particles it is accepted and remains permanently fixed. Otherwise it 
is rejected. In the subsequent desorption step, particles in the adsorbed configuration are 
removed at random until the final density is p2- In the final step, the irreversible adsorption 
process is resumed. New particles can be added to the depleted configuration until the 
jamming limit is reached. As a result of the inherent memory of the RSA process, the 
coverage in this final state, Poo{pii P2) is, in general, different from the usual jamming limit 
for RSA that occurs without the intermediate desorption step, namely 0.547 [§. 

The adsorbed configurations are characterized by a set of distribution functions 
Pn(ri, r2, ■ ■ ■ ; i"n; p), which represent the density of finding n unspecified particles among 
A^ particles defining the system, at the positions ri, r2, ■ ■ ■ , r„ within an area A such that 
p = N/A. Centered on each pre-adsorbed particle is a circular exclusion region of diameter 
a into which it is impossible to insert the center of an additional particle. The probability of 
finding a cavity of diameter 2a centered on ri that is free from the centers of pre-adsorbed 
particles is related to the distribution functions via 

<I>(r?) = l + pjd2fu + ^jd2 rf3/i2/i3P^')(2,3,p) + ^ J d2 d3 d4fufi3fuP^'\2,3,4, p) + 

(1) 

where dl is shorthand for dri etc., and fij = f{rij) = exp{—u{rij)/kT) — 1 is the Mayer 
function. When the interparticle potential, u{rij), is taken for non-overlapping hard spheres 
of diameter a, the Mayer function becomes: 

Equation || applies to any configuration, equilibrium or not, of spherical particles. 

Our strategy is to determine the relationship between the configurations before and after 
the desorption process. This allows us to find an approximate expression for $ which we 
can use to estimate the saturation coverage as a function of pi and p2- 

A. Initial Adsorption 

The kinetics during the irreversible adsorption phase are given by 

'i^K. (3) 

where ka is the adsorption rate constant. 

To solve for the kinetics, one writes down a hierarchy of rate equations for the distribution 



functions [|IO],|Il|. In particular, the pair density evolves according to 

ap(2)(i,2,p) 



dt 

3 



2ka^l,2°) (4) 



where $(1,2°) is the probabihty density function associated with finding an unspecified 
particle at ri and a cavity of diameter 2a free from particle centers at r2. This function may 
also be expressed in terms of Mayer functions and the n-particle density functions: 

$(1, T) = (1 + /i2)[p + / ci3 /23P^'ni, 3, p) + 0{p^)\ (5) 

Substituting ^ in @ to remove the time dependence and noting that p''^-*(l, 3, p) = p^(l + 
/i3)(l + 0(p2)) yields 

gp(^)(l,2,p) _ 'l'(l,2") 
dp $(P) 

_ 2(1 + j\,)[p + p^Jd3 /23(1 + /is) + 0{p')] 
l + p/d2/i2 + 0(p2) 

= 2(l + A2)[p + p'/rf3/i3/23 + 0(p^)] (6) 

Integrating between and p gives the pair density function to third order in p^: 

p(2)(l,2,p) = (1 + f,,)[p' + Ip'JdS /13/23 + 0{p')] (7) 



B. Desorption 

The subsequent desorption process results in the removal of a fraction of the particles 
from an initial configuration at a density pi so that the final density is p2. During this step, 
the density evolves according to 

where k^ is the desorption constant and, if no other relaxation processes (such as surface 
diffusion) are operative during the desorption, the pair density function evolves according 
to 

^^"",f^"^' --2W''(1.2;P..p) (9) 

where we indicate explicitly the dependence of the pair density on the density prior to 
desorption, pi. It follows from these two equations that 

9p(^)(l,2;pi,p) 2 (2).. . . ...s 

= -p^ > (1, 2; pi, p) (10) 

Integrating between the initial (pi) and final densities (P2) yields 

p(2)(l,2;pi,p2) pi 



p(2)(l,2;pi,pi) pI 



(11) 



or 

^(2)(i, 2, pi) = (?(2)(i^ 2, pi, pi) = (?(2)(i, 2, pi, p2) (12) 

where p'-^''(l, 2, p) = p'^g^'^'{l,2, p) defines the radial distribution function. This physically 
intuitive result implies that the structure of an arbitrary configuration is not altered by 
desorption, as long as the particles desorb randomly. Further consequences of this result 
and its connection with quenched annealed systems are discussed elsewhere [^. 

C. Re- Adsorption 

We now wish to relate the properties of configurations during the re-adsorption step to 
those of the configurations in the initial adsorption step. Up to third order in density, the 
differences between the two configurations come from the term p(^^(l,2,p). 

Integrating (|) between p2 and p > p2, 

p(^)(l,2;pi,p2,p)=p(2)(l,2;pi,P2,P2) + (l + /i2)(p' + pD + (l + /i2)(^p'-^P^)/rf3/23/i3 + 0(p^) 

(13) 
Using {\n^ to relate the pre- and post-desorption pair distribution functions, we find 

p(2)(l,2;pi,p2,p) = (1 + /i2)[p' + llp' + pUpi - P2)]Jd3 A3/23 + 0{p')] (14) 

Substituting in (|1|) then yields 

^PuP2,p) = 1+ P J d2f^2 + \p' J d2d3 {1 + f23) (15) 

+ ^[P' + P2(Pl - P2)]Jd2 dS d4 /l2/l3/24/34(l + /23) (16) 

+ ^p3 I d2 d3 dA fuflsMl + /23)(1 + /24)(1 + /34) + Oip') (17) 



where 



$^^^(P)+A32P^(P1-P2), p>P2 (18) 



128/7rV3 9\ 
^-=3^(— -4 (1^1 



is the portion of the 3rd order RSA coefficient due to first order corrections of g^'^^ ||10|| . 
We note that <l>(pi,p2,p) correctly reduces to $^^^(p) in both the limit where pi = p2 (no 
desorption) and P2 = (total desorption). 
During desorption, it is easy to show that 



^p^,p) = <l>^^^(p) + A32p'(pi - p) + 0{p'), Pi > p > P2 (20) 

We plot that $ during adsorption, $i^sA^ desorption $(pi, p), and re- adsorption, $(pi, P2, p), 
in Fig. [VII] . One can also show that 



dp 



_ d<^ 

L.ads dp 



l,des 



+A32PI 



(21) 



where the subscripts ads and des indicate the adsorption and desorption branches, repsec- 
tively and since A32 > 0, we have 



d^ 

dp 



l,ads dp 



2,dcs 



SO the slope of $ during desorption is more negative. One has, furthermore, that 



d^ 
dp 



2,ads dp 



2,des 



-A32P2(2pi - 3p2 



SO 



d^ 
dp 



2,ads dp 



, P2 < -Pl 
2,des O 



(22) 



(23) 



(24) 



d^ 
dp 



2,ads dp 



, P2 > -pl 
2,des O 



(25) 



and thus the change in $ may be more pronounced in either desorption or re- adsorption. 
It is clear that in either case, the available surface function is enhanced compared with the 
simple RSA value (at least to third order in density) since 



$(pi, P2, p) - $^^^(p) = A32plipi - P2) > 0, p > P2 

To obtain quantitative results, we consider the interpolation formula 
<l'(pi,P2,p) = [1 + A32p^(pi -p2)](l-x)^(l + aix + a2x2) 
where 

X = p/Poo 

This form is consistent with the asymptotic kinetics. 



Poo- P ^t 



-1/2 



(26) 



(27) 



(28) 



(29) 



Finally, (|27D correctly reduces to the equation used to determine poo for simple RSA when 
Pl = P2 or when p2 = 0. Expanding the interpolant in powers of p and equating the 
coefficients with those of (0) yields an estimate for Poo(pi5 P2)- As can be seen in Figure |VI1| , 
the saturation coverage is a monotonically increasing function of pi and exhibits a maximum 
versus p2. Actually, the saturation density increases monotonically with the combination 
pKpi — P2), so a maximum occurs for p2 = (2/3)pi. 

In one dimension, the coefficient ^32 takes the value 2/9. Unfortunately, however, the 
analogous interpolation formula, $(p) = [1 + ^32P2(Pi ~ P2)](l — x)^{l + aix + 02^^) does 
not yield physically acceptable estimates of the jamming density, even in the simple RSA 
case. 
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III. SIMULATION 

The adsorption-desorption process was simulated using techniques that have been de- 
scribed previously. Disks are deposited randomly and uniformly in a periodic cell of area 
A, subject to the usual RSA constraints of no overlap and no desorption, until the coverage 
reaches the preset value, pi. In the desorption step, randomly selected particles are removed 
until the coverage falls to p2- The RSA process is then restarted and continued for a total 
of A'triai attempts. In this way we obtain the time-dependent coverage, p(t), where t = Oreii, 
Orel = 7^(T^ /4:A and i is the number of attempts that have been made. This function is then 
averaged over a number, A^av, of independent runs. The coverage in the jamming limit was 
estimated by assuming that at long times the kinetics are described by the usual asymptotic 
kinetics, (^). 

The simulation procedure in one-dimension is somewhat different. Rods of unit length 
are deposited randomly on a line of length L subject to the usual constraints of no overlap 
and no desorption. When the jamming limit coverage is nearly reached, all the remaining 
available gaps, each of which can by now accommodate at most one additional rod, are filled 
by inserting rods randomly and uniformly within them. Partial desorption is then simulated 
by looping through all the adsorbed particles and removing each with a probability pdes- 
This is achieved by generating a uniform random number between zero and one. If this 
number is less than p^es the particle is removed. The depleted configuration then serves 
as the initial configuration for a new irreversible adsorption process which continues until 
jamming. If A^ denotes the number of adsorbed rods, then the coverage is computed from 
p = N/{L — 1). Since the jamming limit can be reached exactly in the one-dimensional 
system, the final results are expected to be more accurate than in two dimensions. 

IV. RESULTS 

In both one and two dimensions, the new jammed state has a higher coverage than 
that of simple RSA, although the effect is smaller than predicted by the theory. In one 
dimension. Figure [VTT| shows that the maximum density (just over 0.78) occurs forp^e^ ~ 0.4 



or p2 ~ 0.6, which is consistent with, although a little larger than, the theoretical prediction 
ofp2 = (2/3)pi. 

In two dimensions (Figure |VII|) , the enhancement is not statistically significant for pi = 
0.4. For pi = 0.53 the effect is clearly evident and the maximum final coverage occurs for 
P2 ~ 0.4 while the theoretical estimate is 0.353. 

That the theory overestimates the enhancement is not totally surprising. The estimate 
of the saturation coverage provided by the interpolation formula for simple RSA (0.553) 
exceeds the actual value (0.547). Moreover, the pair density function is taken only to second 
order in p. 

V. APPLICATION 

An interesting consequence of the findings of this article is that it may be possible to 
achieve a given surface coverage using an intervening desorption process in a shorter time 



than with a continuous adsorption process. We can readily quantify this idea by minimizing 
the time required to reach a given surface coverage in an adsorption-desorption-adsorption 
process by optimizing densities at which the desorption step begins, pi, and ends, p2- Let tf 
denote the time at which the adsorbed density is pf. In a process with a single desorption 
step 



ti 



'^ads ~r t-di 
r-pi 



'breads 
RSAf „\^-l 



/ [^^'\p)]-'dp+ [k,pr'dp+ [<i.(pi,p2,p)]-^rfp 

■^0 J pi J p2 



(30) 
(31) 



We find the minimum time to reach the density pf by by optimizing the density at the 
beginning of the desorption step and that after desorption. As we will demonstrate, the 
amount adsorbed in such a process can exceed that in a continuous adsorption process 
without interruption. The required densities are found from a numerical solution of the 
equations 



dt 



f 



dpi 

dp2 







(32) 

(33) 



This is achieved by iterating the following matrix equation: 



Pi 

P2 



Pi 

P2 



dp'i 



IL. 



IL. 



dpidp2 



dpidp2 
dp'i 



dpi 
dtf_ 

dp2 



The calculations were performed using an interpolation formula. 



$ 



RSA 



(P) 



1 + bix + b2X^ + hx^) 



(34) 



where x = p/poo and 6i, 62, and poo are taken from reference ||T0[- Fig. |VI1] shows pi and p2 



which minimize the time required to produce an adsorbed configuration with a final density 
of Pf. For a given value of kd, the two curves converge at p/o. If the desired final density is 
larger than this value, it will be more efficient i.e., require less time, with a two-step process 
than with a single step process. If the final density is less than p/o, the single step process 
is more efficient. As k^ increases, pjQ decreases since the time required to remove a given 
fraction of particles is a decreasing function of k^. Thus at large desorption rates less time 
is spent in the desorption phase. We observe that pi appears to be almost independent of 
kd (for Pf > p/o), while p2 curves converge at large values of kd. The two-step process is 
most efficient in the limit fc^ — > 00. That the two step process can be more efficient than a 
one-step process is a result of the inherent memory of the adsorbed configurations. Without 
this memory effect, the two-step process is always less efficient. 

VI. CONCLUSION 

We have developed a theoretical approach for modeling a three-step surface filling process 
composed of two irreversible adsorption steps separated by a desorption step. This approach 
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considers the effects of surface density and coverage history on the adsorbate particle dis- 
tribution functions and predicts an enhancement of the saturation surface density that is 
maximized when one third of the molecules are removed during the intermediate desorption 
step. 

One experimental realization of this process is partially reversible macromolecular ad- 
sorption where only a certain fraction of molecules desorb during a buffer rinse. Partial 
reversibility may be due to adsorbent heterogeneity, formation of stable clusters of molecules 
on the surface, or post-adsorption transitions in conformation or orientation. Another ex- 
perimental realization of this proces is the elution of otherwise irreversibly adsorbed macro- 
molecules with surfactant or detergent eluting agents [^. The results presented here could 



have important implications for these and other multi-step filling processes when a high 
saturation density and/or a rapid surface filling is required. 
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FIGURES 

FIG. 1. The available surface function plotted as a function of density at various fixed densities 
prior to and following desorption. The bottom line is properly placed and each subsequent line is 
shown displaced vertically by 0.1 units for clarity. 

FIG. 2. The jamming limit density, as determined from the interpolation formula ([27|), as 
a function of (a) the density, pi at which desorption begins and (b) the density, p2, at which 
desorption ends. 

FIG. 3. Enhancement of the jamming limit coverage, pf, in one dimension resulting from 
a random desorption from a jamming RSA configuration of hard rods on a line. Each particle 
is desorbed from this configuration with probability pdes before re-adsorption to a new jamming 
limit. 

FIG. 4. Enhancement of the jamming limit coverage, pj, in two dimensions resulting from a 
random desorption at a density pi to a density p2 as determined in a numerical simulation of the 
process. +: pi = 0.53, *: pi = 0.5, diamonds: pi = 0.4. 

FIG. 5. Kinetics of an adsorption-desorption-adsorption process, (a) For a given desorption 
rate, k^i, the values of the densities pi and p2 that minimize the time required to reach the final 
coverage p/.(b) Comparison of the optimized two step process with k^ = 1 and optimized for 
Pf = 0.5 with an uninterrupted adsorption process. Note that the final coverage is attained more 
rapidly using the two-step process. 
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TABLES 

TABLE L Enhancement of the Saturation Coverage. Particles are desorbed randomly from an 
RSA configuration at an initial density pi until the density reaches p2- A new RSA process then 
begins on the depleted configuration. The density at saturation is then estimated by extrapolating 



the coverage versus time curves using the asymptotic power law, (29). Each condition used 800 
independent averages for a system with a relative area of the particle to simulation cell of 0.001. 
For 9i = 0.53 a total of 10^ attempts were made to add new particles, while for 6i = 0.50 and 
9i = 0.4 we used 5 x 10® attempts. 
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